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ACCEPTABLE TOLERANCES FOR MATCHING ICING SIMILARITY PARAMETERS 

IN SCALING APPLICATIONS 

David N. Anderson* 

Ohio Aerospace Institute 
Brook Park, Ohio 44142 


Abstract 

This paper reviews past work and presents new data to 
evaluate how changes in similarity parameters affect ice 
shapes and how closely scale values of the parameters 
should match reference values. Experimental ice 
shapes presented are from tests by various researchers 
in the NASA Glenn Icing Research Tunnel. The 
parameters reviewed are the modified inertia parameter 
(which determines the stagnation collection efficiency), 
accumulation parameter, freezing fraction, Reynolds 
number and Weber number. 

It was demonstrated that a good match of scale and 
reference ice shapes could sometimes be achieved even 
when values of the modified inertia parameter did not 
match precisely. Consequently, there can be some 
flexibility in setting scale droplet size, which is the test 
condition determined from the modified inertia 
parameter. A recommended guideline is that the 
modified inertia parameter be chosen so that the scale 
stagnation collection efficiency is within 10% of the 
reference value. The scale accumulation parameter and 
freezing fraction should also be within 10% of their 
reference values. The Weber number based on droplet 
size and water properties appears to be a more 
important scaling parameter than one based on model 
size and air properties. Scale values of both the 
Reynolds and Weber numbers need to be in the range of 
60 to 160% of the corresponding reference values. The 
effects of variations in other similarity parameters have 
yet to be established. 

Nomenclature 

A c Accumulation parameter, dimensionless 

b Relative heat factor, dimensionless 

c Airfoil chord, cm 

c P ,ws Specific heat of water at the surface 
temperature, cal/g K 

d Cylinder diameter or twice the leading-edge 

radius of airfoil, cm 


*Senior Research Associate, Member A1AA 


h c Convective heat-transfer coefficient, 

cal/sm 2 K 

h G Gas-phase mass-transfer coefficient, g/sm 2 

K Inertia parameter, dimensionless 

K 0 Modified inertia parameter, dimensionless 

LIVC Cloud liquid-water content, g/m 3 

M Mach number, dimensionless 

MVD Water droplet median volume diameter, pm 

n Freezing fraction, dimensionless 

p st Static pressure, Nt/m 2 

r Recovery factor, dimensionless 

Re Reynolds number of model, dimensionless 

Res Reynolds number of water droplet, 

dimensionless 

s Distance from stagnation line, cm 

tf Freezing temperature, °C 

t s Surface temperature, °C 

t s , Static temperature, °C 

t lo , Total temperature, °C 

V Air velocity, m/s 

We Weber number based on droplet size and water 
properties, dimensionless 

We c Weber number based on model size and air 
properties, dimensionless 

P Local collection efficiency, dimensionless 

po Collection efficiency at stagnation line, 

dimensionless 

6 Air energy transfer parameter, °C 

A Droplet range, m 

Astons Droplet range if Stokes Law applies, m 

Af Latent heat of freezing, cal/g 

A v Latent heat of condensation, cal/g 

p Air viscosity, g/m s 

ju w Liquid water viscosity, g/m s 

p Air density, g/m 3 

Pi Ice density, g/m 3 

p w Liquid water density, g/m 3 

a Surface tension of water over air, dyne/cm 

T Accretion time, min 

(f) Droplet energy transfer parameter, °C 

Subscripts 

R Reference 

S Scale 
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Introduction 

Although a number of similarity parameters have been 
identified as possibly being important in determining 
ice shape, there have been only a few studies that have 
attempted to evaluate how each parameter affects the 
ice shape or how closely scale and reference values of 
the similarity parameters must agree. This paper 
reviews both published and unpublished ice-shape data 
in an attempt to gain further insight. 

By the 1970’s it was well recognized that there were at 
least 6 basic similarity parameters that needed to be 
considered in performing scaling tests. They were the 
modified inertia parameter, K 0 , defined by Langmuir 
and Blodgett 1 , the accumulation parameter, A c , the 
freezing fraction, n, the water-energy transfer 
parameter, cj), the air-energy transfer parameter, #and 
the relative heat factor, b. The parameters n, </>, and 6 
are from the analysis of Messinger, 2 and b was defined 
by Tribus. 3 Different scaling methods match scale 
values of various parameters to their respective 
reference values to obtain a set of simultaneous 
equations that are solved for the required scale test 
conditions. Charpin and Fasso 4 were the first to 
incorporate most of these parameters into a scaling 
method for sea-level testing facilities. By the mid 
1980’s Ruff 5 established that the most faithful scaling 
resulted when five of these parameters (ignoring b) 
were matched in a facility with control of test-section 
pressure. Soon after, Bartlett 6 7 8 performed a series of 
analytical and experimental studies to try to determine 
how important these parameters were to the ice shape. 
Could any of them be ignored? How closely did each 
have to be matched to insure an acceptable scaled ice 
shape? 

There are several important reasons to do studies of the 
kind Bartlett did. First, in planning scale tests, it is not 
always possible to achieve precisely the scaled 
conditions required to match all of the scaling 
parameters. Therefore, some flexibility in choosing test 
conditions is helpful. For example, if the test facility 
cannot produce the droplet size needed to exactly match 
modified inertia parameters, it is useful to know if other 
values inside the tunnel operating map might provide 
ice shapes with acceptable agreement. Second, the test 
conditions are only known within some uncertainty 
band; thus, it is important to know what levels of 
uncertainties are acceptable when doing scaled tests. 
Finally, identifying the similarity parameters with the 
greatest effect on ice shape can be helpful in 
understanding the physics of icing. 

The questions that Bartlett began to consider have not 
yet been adequately answered. The present approach to 


scaling is still to match each parameter precisely and 
determine the scaled test conditions accordingly. There 
is now, however, additional data which will be 
reviewed in this paper to add further insight into the 
relevance and acceptable variation of several similarity 
parameters. 


Similarity Parameters 

The similarity parameters used in icing scaling have 
been identified by examining the physics of icing. 
Similarity of each of several processes needs to be 
maintained between the scale and the reference 
situations. 

The first of these processes is the flowfield. It can be 
simulated by matching scale and reference values of Re 
and M, by using a model that is dimensionally similar 
to the full-scale (reference) article and by setting the 
scale angle of attack equal to the reference. For icing 
encounters, the speeds involved are typically low 
enough that compressibility is not an issue, and M is 
neglected. Although Re has usually been ignored as 
well, Bilanin 6 advocated including it in icing scaling 
analyses, and some recent studies 10 have used Re = 
Vdp/jU as an optional scaling parameter. Here d is the 
diameter for cylindrical models or twice the leading- 
edge radius for airfoils. In this study, d was also used 
as the geometric length scale for the parameters K, A ,,, 
h c , and We c . For the NACA 0012, a leading-edge 
radius of .0158c was used; for the business jet, .013c, 
where c is the model chord. 

The next process to be represented is the droplet 
motion. Similarity of droplet trajectories can be 
obtained by matching the modified inertia parameter, 
K„, of Langmuir and Blodgett: 1 


K 



A 

^ Stokes 



( 1 ) 


AJAstokes is the droplet range parameter, defined as the 
ratio of actual droplet range to that if Stokes law for 
solid-sphere drag applied. The range parameter is a 
function only of the droplet Reynolds number. Res 
Langmuir and Blodgett presented tabulations of this 
function. K is the inertia parameter, 


Pw MVP 2 V 
18 dp 


( 2 ) 


The water-droplet collection efficiency at the stagnation 
line is directly dependent on the modified inertia 
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parameter. Langmuir and Blodgett gave this 
relationship: 

_ 1.40(W 0 -l/8f 84 

“ l + 1.40(^ 0 -l/8) 0 84 (3) 

Eq. (3) was used to calculate all the j5 n ' s for this report. 
Although this relationship gives only the value of 
collection efficiency at the stagnation line, if two 
similar airfoils have the same K 0 , the collection 
efficiency will also match everywhere on the two 
airfoils. 


rime ice, water freezes on impact so that there is no 
liquid-water film; therefore, surface dynamics are only 
relevant to glaze-ice accretions. The physical processes 
on the liquid surface are not yet well understood, but 
similarity parameters related to surface effects can be 
identified that may be important. These include the 
capillary number, N cap = jU K VI a, included by Kind" in a 
recent study, as well as the Weber numbers, We = 
V 2 MVD pjo and We c =V 2 dpto. Kind suggested that 
a Weber number based on the water-film thickness 
might also be considered. In addition, the free-stream 
Reynolds number, Re, probably plays a role in water- 
film dynamics. 


The third consideration is total ice accretion. The 
quantity of ice is represented by the accumulation 
parameter, A c : 


LWCVt 
d Pi 


(4) 


A match of scale and reference A c insures that the same 
quantity of ice, relative to airfoil size, is accreted for 
both. This matching assumes that J3 0 for scale and 
reference situations also matches. 

Next, for glaze ice, similarity in energy balance is 
required. Messinger’s 2 surface energy balance can be 
written in the form, 


c p,WS 




(5) 


where n is the freezing fraction; <]). the water energy 
transfer parameter; 6, the air energy transfer parameter; 
and b, the relative heat factor, introduced by Tribus, et. 
al. 3 These parameters are defined as 


4> = tf -*st 


V 2 


(6) 


Finally, the non-dimensional water-film thickness itself 
may well be a similarity parameter of importance in 
describing water-film dynamics. Feo and Urdiales 12 
gave an experimental correlation for water-film height 
for heavy-rain conditions, and Feo 13 also recently 
measured the film thickness for near-Appendix-C cloud 
conditions in a warm-air tunnel. No other experimental 
studies applicable to icing are known to have been 
published. A limited unpublished study by Anderson 14 
showed good icing scaling results when Feo’s latest 
expression for film thickness normalized to the model 
size was matched between scale and reference 
conditions. While encouraging, these results are still 
preliminary and will not be included here. 

In this paper, the sensitivity of ice shapes to variations 
in K 0 , A,, n. Re, We and We, will be discussed. Other 
parameters will not be considered because there is 
insufficient data to evaluate their effect. K„ and n are 
determined at the stagnation line of the model only. It 
is assumed that values there will be proportional to 
values anywhere on the airfoil. The effects of static 
pressure, p sl , on both the similarity parameters and ice 
shape will also be reviewed to help assess the 
importance of some of the parameters for icing. The 
evaluations will be made for glaze ice. The ice shapes 
used for illustration were obtained in the NASA Glenn 
Icing Research Tunnel by various researchers. 


9 = 


\s *st 


rV l 

2c 


+ - 


P WW Pw 


k (7) 


p J 


b = 


LWC V Pq c p 


p,ws 


( 8 ) 


The last processes to be considered occur on the surface 
of the model where water collects before freezing. For 


Effects of Parameters on Ice Shape 
Modified Inertia Parameter. K „ 

An indication of the possible effect of K 0 on ice shapes 
can be deduced from its effect on collection efficiency. 
Figure 1 shows the relationship between J3 0 and K„ 
given by eq. (3). Note that for K 0 greater than 1, f} 0 
increases less than 10% for a K 0 increase of 25%. 
Between 2 and 3, K 0 has to vary by more than 50% for 
P 0 to change 10%. When K 0 has reached a value of 3, 
the collection efficiency is already nearly 80%, and K 0 
has to triple to a value of 9 for /3 0 to increase from 80 to 
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Modified Inertia Parameter, Kg 

Figure 1 . Relationship between K 0 and the Langmuir 
and Blodgett 1 /?,,. 

90%. Bartlett 8 assumed that a ±10% variation in ice 
shape could typically be expected when tests were 
repeated. Therefore, a 10% variation in /?„ is probably 
within run-to-run variations, and it is apparent that over 
much of the range of K 0 shown in figure 1, a perfect 
match of this similarity parameter is unnecessary to 
insure scaling results within repeatability expectations. 

K n and collection efficiency are dependent on velocity, 
droplet size and model size. Of most interest to 
practical scaling applications is how much flexibility in 
the choice of scale MVD results from relaxing the 
requirement that K os = K nR . While it is always 
preferable to test with scale droplet sizes as close to the 
exact solution as possible, in some situations, the 
scaling equations may yield scale MLD’s that are 
outside the tunnel cloud capability. Estimates of 
permissible variations in droplet size can be made from 
plots of P 0 as a function of MVD. Figure 2 shows two 
sets of such plots. Each curve represents either a 
reference or scale condition for an NACA 0012 airfoil 
of a specific size. A reference velocity of 67 m/s was 
chosen, and the scale conditions were determined using 
the Ruff method, described in reference 10. The scale 
velocity was the average of the velocities that resulted 
from applying constant- We and constant-Re constraints. 


and solid circles, respectively. The stagnation-line 
collection efficiency for these drop sizes was .80. For 
the example shown, it can be seen that a scale droplet 
size as large as 25 pm would only increase the scale 
stagnation-point collection efficiency by 10%. 

Larger models allow less flexibility in choosing the 
scale droplet size. For example, figure 2(b) gives the 



MI D , pm 


(a) Scaling from 53.3- to 26.7-cm Chord. 



MVD. pm 


Figure 2(a) shows curves for 53- and 27-cm-chord 
airfoils. For a reference MVD of 30 pm, a scale MVD 
of 15.6 pm was found by matching scale and reference 
K 0 . The reference and scale MVD’s are shown as open 


(b) Scaling from 183- to 91-cm Chord. 

Figure 2. Examples of Collection Efficiency for 14- 
Size Scaling. NACA 0012 Airfoils at 0°AOA. 
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same scaling situation as figure 2(a), but for reference 
and scale model chords of 183 and 91 cm. Although 
the scale MVD found from K os = K 0 ,r is still 15.6 pm, 
the collection efficiency for both models is now only 
.59. Restricting the scale J3 0 to be no more than 10% 
higher than the reference requires that the scale drop 
size be no larger than 19 pm, instead of the 25 |im 
found for the conditions of figure 2(a). 

If alternate drop sizes are contemplated, the full 
collection-efficiency curve over the region of expected 
ice accretion for the model at the off-scale conditions 
should be compared with that for the reference case to 
insure that differences are small. It is also important to 
remember that measured drop sizes cannot be precisely 
defined. In the 1RT, for example, the drop size 
uncertainty has been estimated to be ±11%. 15 Thus, the 
difference between the alternative drop size and the 
exact-solution value must be greater than 11% to be 
significant. 

This flexibility in choosing scale drop size can be 
confirmed with experimental results. Ice shapes 
obtained in the IRT by various researchers are 
compared in figures 3 and 4. The ice-shape data were 
reported by Wright. 16 All were obtained on 53 -cm- 
chord NACA 0012 airfoils at 4° angle of attack. 
Combinations of test conditions were selected for each 
of the comparisons such that constant A c , n, b, (f), and 0 
were maintained, but Kg s were different. 
Coincidentally, Re and We c also matched. The test 


conditions for each test and the corresponding 
similarity parameters are shown below the data plots. 

In figure 3(a), ice shapes produced with drop sizes of 
15 and 40 pm are compared. The corresponding Kg s 
were 1.63 and 2.56. The jig were calculated to be .66 
and .75, respectively. In addition to the K 0 differences, 
the We for the two cases were 2.44 and 3.25xl0 3 . The 
ice shapes differ mainly in the size and extent of the 
upper horn, with the higher K n leading to a slightly 
larger horn. The differences were little more than 
typical run-to-run variations. 

Figure 3(b) compares the clean-airfoil collection- 
efficiency curves for the conditions of the two icing 
tests of 3(a). The /? s were calculated using LEWICE 
2.0. 17 The largest differences between ji curves were 
seen on the lower surface. In spite of these differences, 
the lower-surface ice shown in figure 3(a) agreed well 
for the two tests, presumably because the absolute value 
of ji was low where the curves deviated the most. The 
ji curves showed only minor differences between the 
accretion limits for the two tests. 

Figure 4(a) compares shapes formed withMFD’s of 25 
and 40 pm. These drop sizes corresponded with K 0 = 
3.62 and 7.54. The J3 n were .80 and .88, respectively. 
Again, all other similarity parameters were matched 
except We. The ice shapes were in close agreement 
including both the upper horn and a large feather 
structure just aft of the upper horn. It can be seen from 
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Run 
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v. 

MVD, 
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T, 

Kg 

Po 

A 

n 

U 

<P, 

e. 

Re, 

We. 

We c , 

°C 

°c 

m/s 

pm 

g/m 3 

sec 


0 

°C 

°c 

10 4 

10 3 

10 3 

314 

-10.8 

-5.6 

102.8 

15.0 

0.60 

384 

1.63 

0.66 

1.53 

0.53 

0.34 

9.6 

11.2 

13.0 

2.44 

3.40 

303 

-10.8 

-5.6 

102.8 

20.0 

0.55 

420 

2.56 

0.75 

1.54 

0.52 

0.36 

9.6 

11.2 

13.0 

3.25 

3.40 


Figure 3. Effect of K n on Ice Shape. 53. 3-cm-Chord NACA 0012 Airfoil at 4°AOA. Ice-Shape Data from Addy, 
April 1997, reported in CD’s accompanying Wright. 16 Collection Efficiencies Calculated Using LEWICE. 17 
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x/c 

MVD = 25 urn. Run 414 

(a) Ice Shapes 


Run 


t/oh 

v. 

MVD, 

LWC, 

% 

K 0 

°c 

°c 

m/s 

pm 

g/m 3 

sec 

414 

-10.8 

-5.6 

102.8 

25.0 

0.55 

420 

3.62 

413 

-11.2 

-5.9 

102.8 

40.0 

0.50 

462 

7.54 



(b) Collection Efficiency 


Pn 

A c 

n 

b 

</>, 

°C 

e, 

°c 

Re, 

10 4 

We, 

10 3 

We a 

10 3 

0.80 

1.54 

0.49 

0.38 

9.6 

11.2 

13.0 

4.06 

3.40 

0.88 

1.54 

0.51 

0.38 

9.9 

11.7 

13.0 

6.50 

3.41 


Figure 4. Effect of K 0 on Ice Shape. 53. 3-cm-Chord NACA 0012 Airfoil at 4°AO A. Ice-Shape Data from Bidwell 
and Van Zante, January 1998, reported in CD’s accompanying Wright. 16 Collection Efficiencies Calculated Using 
LEWICE. 17 


figure 4(b) that the lower-surface collection-efficiency 
for the two runs did not match even as well as those for 
the cases given in figure 3, but again these differences 
did not appear to have much effect on the lower-surface 
ice actually accreted. 

These experimental comparisons show that scale and 
reference values of /?« can differ by 10 - 1 1% without a 
significant effect on ice shape. As a guideline, when 
the scale droplet size determined by setting K os = K ll R 
is outside the capability of a test facility, alternate 
values of MVD can be substituted, providing for the 
new MVD is within 10% of the reference value and p 
curves for the scale and reference airfoils do not deviate 
significantly. 

Accumulation Parameter, A 

The accumulation parameter directly determines the 
quantity of ice produced. In scaling calculations, the 
matching of scale and reference A c (see eq. (4)) allows 
the determination of the scale spray time. In principle, 
time can be set easily in test facilities, so that there is no 
need for the scale value of this parameter to deviate 
from the reference. However, for short spray times, 
spray-system stabilization may be a significant portion 
of the total spray for some facilities. The IRT now has 
a rapid-start spray system that avoids this problem. The 
uncertainty in A c has been estimated to be about ±12% 15 
in the IRT, due to uncertainties in the test conditions 


other than time. Thus, the uncertainty in the total 
quantity of ice accreted will be no better than ±12%, 
and to avoid misleading results, shielding of the model 
during spray stabilization or corrections to the accretion 
time may need to be considered for some facilities. 

If the scale stagnation collection efficiency is different 
from the reference, instead of matching A c , the product 
A c p 0 should be matched as closely as possible between 
reference and scale. 

Freezing Fraction, n 

The freezing fraction is primarily a function of 
temperature, airspeed and L WC, but it is the scale L WC 
that is usually found by equating scale and reference 
freezing fraction. For most situations, there is no 
practical advantage in allowing the scale n to deviate 
from the reference, for if the scale value of L WC found 
initially from solving the scaling equations is not 
suitable, an alternate value can be selected by applying 
the Olsen method. 18 This method requires simply that 
the scale temperature be adjusted for the new LWC so 
that n s = n R is maintained.. However, due to test- 
condition uncertainties, n is only known to within about 
±10% 15 , so it is useful to know what kinds of variations 
in ice shape can be expected. 

The importance of the freezing fraction in determining 
glaze ice shapes has long been recognized. In 1986, 
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Olsen, et al. 19 showed the effect of systematic changes 
in temperature on ice shape with all other test 
conditions fixed. When freezing fractions are 
calculated from Olsen’s reported test conditions, it can 
be shown that ice shapes changed little when n varied 
by less than 10%. In reference 18, tests were reported 
in which LWC was varied and ice shapes were recorded 
on cylinders mounted vertically in the IRT. When the 
temperature was adjusted to maintain constant /?, ice 
shapes showed only minor random variations. When 
temperature was fixed so that n varied with LWC, the 
ice shapes changed more significantly. 

Some of the results from reference 18 will be repeated 
here along with other unpublished data from that study. 
Figure 5 gives ice shapes resulting when the freezing 
fraction was varied with all test conditions except LWC 
fixed. Note from the table accompanying figure 5 that 
all the similarity parameters except n and b were also 
constant. For the test conditions in this example, an 
increase in LWC from .8 to .9 g/m 3 reduced the freezing 


fraction from .38 to .35. Figure 5(a) shows how the 
resulting glaze-ice horns spread apart slightly, shifting 
aft on the model as n decreased. The change in horn 
angle for this modest change in n was no greater than 
that typically seen when tests are repeated at the same 
test conditions. A further increase in LWC from .9 to 
1.0 g/m 3 reduced n from .35 to .33 and resulted in a 
slight additional spreading of the horns. The overall 
change in n was about 13%, which produced a 
noticeable, but probably not significant, effect on ice 
shape. 

Figure 5(b) gives the results of larger changes in n. 
LWC’s were .8, 1.1 and 1.3 g/m 3 with corresponding 
freezing fractions of .38, .31 and .29. Consistent with 
the trend seen in figure 5(a), the horns spread out as the 
freezing fraction decreased. The reduction of n from 
.38 to .29, representing a change of 24%, produced 
significant changes in horn position for these fairly 
large accretions. For small accumulation parameters 
(short accretion times), however, differences over the 



x/c 


1-13-98/18; /? = 0.38 

1-13-96/19; ;? = 0.35 

1-13-93/20; /? = 0.33 



x/c 


1-13-93/18; n = 0.38 

1-13-93/21;;? = 0.31 

1-7-93/13; n = 0.29 


(a) Liquid-Water Contents of 0.8, 0.9 and 1.0 g/m 3 . 
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Figure 5. Effect of Freezing Fraction on lee Shape. 5.1-cm-Diameter Vertical Cylinder Tested in the NASA Glenn 
IRT. Published ls and Unpublished Ice-Shape Data from Tests by Anderson. 
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range of freezing fraction tested in figure 5(b) would 
probably not be substantial. 

The results in figures 5(a) and (b) demonstrate how 
sensitive ice-shape is to freezing fraction. For scaled 
tests for which the objective is to simulate a reference 
ice shape, the ±10% uncertainty in n leaves no 
flexibility to permit n s to deviate from n R . 

Reynolds and Weber Numbers, Re, We and We r 

Any of these parameters can be matched to their 
reference values to determine the scale velocity. While 
they have not usually been included in scaling studies, 
Anderson and Ruff 10 showed that scaling methods 
incorporating Re and We were more successful than 
methods that didn’t. Reference 15 found the best 
scaling when V s was determined by a compromise 


between constant Re and constant We. This 
compromise put both parameters in the range of 
approximately 60 to 160% of the reference values. 

Figures 3 and 4 above compared ice shapes from pairs 
of tests for which Re and We c matched. We for run 303 
in figure 3 was 133% of that for run 314. For the tests 
of figure 4 We for run 413 was 160% of that for run 
414. In both cases, the ice shapes matched well in spite 
of these We differences. These results and those of 
reference 15 suggest that, if We matters at all, an 
approximate match between reference and scale values 
of We is sufficient for good matches of ice shape. 

One of the properties included in the definition of We is 
the surface tension. To demonstrate the importance of 
surface tension, and possibly We, on ice shape, Bilanin 
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Conditions Adjusted to 
Approximate We. 


(c) Model Size Varied and Test 
Conditions Adjusted to Match We 
and Re. 
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Figure 6. Effect of Surfactant on lee Shape. Vertical Cylinders Tested in the NASA Glenn IRT. Published 20 and 
Unpublished Ice-Shape Data from Tests by Bilanin and Anderson. 
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and Anderson 20 performed a series of tests in the IRT in 
August 1993 with surfactant added to the spray water. 
Figure 6 gives three sets of ice-shape comparisons 
resulting from these tests. 

Figure 6(a) compares ice shapes obtained from two 
tests with the same model size and test conditions, 
except that for the first test the spray was with water 
only while the second used a water-surfactant mixture. 
The two tests produced matching values of all similarity 
parameters except for the two Weber numbers. We and 
We c for the surfactant-addition test were more than 
double those for the water-only test. The dramatic 
difference between the two ice shapes shows that 
simply matching K 0 , A c , n, b, <f), 6 and Re is not 
sufficient, and either We or We, also has a significant 
effect on ice shape. 

In figure 6(b) another pair of ice shapes are compared 
with the same model size. Again, one shape was 
obtained with water only and one with surfactant added 
to the spray. For these two tests, however, test 
conditions differed in such a way that We and We c for 
the surfactant-addition test were no more than 124% of 



x/c 

3-13-98 Run 736.31; c = 91.4 cm 

3-6-98 Run 718.31; c = 45.7 cm 


those for the water-only test. The scale Re for the 
surfactant-addition test was about 70% of that for the 
water-only test, and all other similarity parameters 
either matched or nearly matched. While the two ice 
shapes of figure 6(b) were not in perfect agreement, 
they were similar in most characteristics. 

Finally, in figure 6(c), both test conditions and size 
were altered when surfactant was added. Both We and 
Re for the surfactant-addition test were within 10% of 
the water-only test, while We c for the two tests were 
within 21%. Most of the other similarity parameters 
except K q , ^and # agreed within 10% for the two tests, 
as well. Ice shapes were again in close agreement. 

The comparison of the results of figure 6(a) with either 
(b) or (c) provides a powerful argument for including 
We or We c as a scaling parameter. It is clear that while 
scale Re, We and We c may not have to match their 
respective reference values closely, differences in the 
Weber numbers of a factor of 2 are not acceptable. 

These conclusions can be further confirmed by looking 
at the results of a series of 1998 scaling tests performed 
by Chen 21 without surfactant addition. Figure 7 gives 



x/c 

3-13-98 Run 736.31; c= 91.4 cm 

3-5-98 Run 718.21; c = 45.7 cm 


(a) Constant- Velocity Scaling; Re s /Re R = 0.49; (b) Scaling with V S /V R = 1.15; Re s /Re R = 0.58; 

We s / We R = 0.61. We s /We R = 0.82. 
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Figure 7. Scaling to Vi Size. Business Jet Airfoil at 0° AOA Tested in the NASA Glenn IRT. Ice-Shape Data from 
Tests by Chen. 21 
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Chen’s results of scaling from a reference model of 
91 -cm chord to a scale with 46-cm chord. In figure 7(a), 
ice shapes are compared for the situation in which the 
scale velocity (run 718.31) nearly matched the 
reference (run 736.31). Other conditions were such that 
the similarity parameters K 0 , A c , n, <p, and 6 also 
matched. However, the scale Re was about 50% of the 
reference value, the scale We was about 60% and the 
scale We c was about 48%. While the quantity of ice 
obtained in the scale test was about right, the shape did 
not agree with the reference. 

For the shapes shown in figure 7(b), the reference test 
was again 736.31, as in figure 7(a). Scale conditions 
now gave Re, We and We c that were 58, 82 and 67%, 
respectively, of their reference values while K 0 , A,, n, b, 
and # matched. While the scale shape was not a perfect 
representation of the reference, its simulation was 
greatly improved over that of figure 7(a). 

For the test conditions reported in figure 7(c), the scale 
Reynolds number was now 67% of the reference value, 
and We and We c were 104 and 88% of the reference 



- 0.05 0.00 0.05 0.10 


x/c 

3-11-98 Run 736.11; c = 91.4 cm 

3-6-98 Run 718.31; c = 45.7 cm 

(c) Constant-lfe Scaling; V S /V R = 1.33; Re s /Re R = 
0.66; We s /We R = 1.03. 


values, respectively. Simulation of the reference shape 
was quite good except for the upper horn position. 
Scale conditions for these tests again resulted in a 
match of reference values of K 0 , A c , n, b, and 6. 

Finally, for figure 7(d), the scale Re, We and We c were 
73, 124 and 109% of the reference values, respectively, 
and the scale parameters K 0 , A c , n, and # matched the 
respective reference values. The scale ice shape for this 
test was an excellent simulation of the reference shape. 
Without changing surface tension, it is not possible to 
match both Re and We simultaneously. However, the 
scale values of Re and We for the tests of figure 7(d) 
represent a compromise between perfect matching of 
either parameter with their reference values. 

Chen’s data, those of figures 3 and 6 and the results in 
reference 15 consistently indicated that the best scaling 
can be obtained when both the scale Re and scale We 
are between approximately 60 and 160% of their 
respective reference values. 
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Figure 7. (concluded). 
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Effect of Pressure on Similarity Parameters 

A knowledge of the effect of pressure on both ice 
shapes and similarity parameters can help to determine 
which similarity parameters have the most influence on 
ice shape. For example, if a change in pressure 
produces a large effect on a particular similarity 
parameter, but little effect on ice shape, it can be 
concluded that ice shapes are insensitive to that 
parameter. 

Studies to evaluate the effect of pressure on ice shapes 
have been reported by Bartlett 6,7 over the range of 3.0 to 
9.8xl0 4 Nt/m 2 and by Oleskiw, et al.~ 2 for pressures of 
4.6 to lO.lxlO 4 Nt/m 2 . Within the typical repeatability 
of test facilities there was no evidence from either study 
that pressure change affected the ice shape. 

Table I shows how pressure affects calculated values of 
each of the scaling parameters. An NACA 0012 with a 
chord of 53.3 cm was used in the computations. The 
greatest effect of pressure is on K„, 0, b. Re and We c , 
while J3 n , A c , </>, n, and We are not affected, or the effect 
is slight. Even for a chord of 183 cm, the effect of 
pressure on fi n , over the range of pressure shown and 
for the velocity given, would be less than 10%. For 
higher velocities, the variation of p n with pressure will 
be even less. Thus, the lack of an effect of pressure on 
ice shape is consistent with the small effect of pressure 
on modified inertia parameter and stagnation collection 
efficiency. Although both 9 and b change with 
pressure, the ratio 61b is nearly constant; thus, no 
conclusions can be reached about the importance of 
either of these parameters on ice shape. 

Re is nearly directly proportional to pressure due to the 
effect of pressure on air density. This strong pressure 

Table I. Effect of Pressure on Similarity Parameters 
53.3-cm-Chord NACA 0012 Airfoil 

t s „ -6.7°C; t tot , -4.4°C; V, 67 m/s; MVD, 30pm; 

LWC, 1 g/m 3 ; time, 7.3 min. 


p st , 10 4 Nt/m 2 

9.72 

8.43 

6.97 

5.72 

4.65 

K„ 

3.70 

3.90 

4.18 

4.48 

4.80 

Po, % 

80.3 

81.0 

81.9 

82.8 

83.6 

A c 

1.91 

1.91 

1.91 

1.91 

1.91 

n 

0.27 

0.27 

0.27 

0.27 

0.27 

b 

0.57 

0.61 

0.68 

0.76 

0.85 


6.1 

6.1 

6.1 

6.1 

6.1 

0, °C 

8.5 

9.1 

10.1 

11.4 

13.0 

Re, 10 4 

8.54 

7.41 

6.13 

5.03 

4.09 

We, 10 3 

2.07 

2.07 

2.07 

2.07 

2.07 

We c , 10 3 

1.47 

1.28 

1.06 

0.87 

0.71 


effect suggests that ice shapes are fairly insensitive to 
Re, and changes by a factor of 2 or more in this 
parameter appear to be tolerable. Further evaluation of 
Re effects are needed. 

Earlier discussion showed that when the Weber number 
changed by a factor of 2, as We c does for the range of 
pressures given in Table I, the shape of the ice changed 
dramatically. Thus, for ice shapes to be independent of 
pressure, We c cannot be the correct form of Weber 
number. Either We, which is unaffected by pressure 
change, or any Weber number using water properties, 
appears to be the proper expression of this similarity 
parameter. 

Summary 

This study of how ice shapes are affected by variations 
in the scaling similarity parameters has resulted in the 
following tentative guidelines on acceptable deviations 
of scale parameters from the corresponding reference 
values: 

1. The scale K 0 does not need to match the reference 
value closely, but the permissible deviation 
depends on the reference value. It is recommended 
that as a guideline, the scale J3 n should be within 
10% of the reference value, where fi„ is related to 
K 0 according to the expression given in eq (3). 
Furthermore, if a scale MVD is selected that does 
not satisfy K os = K nR , the scale and reference 
collection efficiencies along both the reference and 
scale airfoil surfaces should be computed and 
compared to insure that differences are minor. 

2. Uncertainties in A c should not be permitted to 
exceed about ±10% to insure a reasonable match of 
ice shapes. This similarity parameter depends on 
the spray time. When spray times are short, they 
may need to be corrected to account for spray-bar 
system stabilization periods. When the scale p 0 
does not match the reference value, the product 
A c p o should be matched. 

The following comments pertain only to glaze 
conditions. For rime conditions, none of the parameters 
n. Re and We or We c have an influence on ice shape. 

3. The scale n should be held to within 10% of the 
reference value. This tolerance is the same as 
estimated uncertainties in the freezing fraction. 
Thus, the goal should be to try to match this 
parameter to its reference value. 

4. Acceptable scale Re and We fall in the range of 60 
to 160% of the respective reference values. 
Evaluation of pressure effects on ice shape 
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suggested that even greater tolerances on Re may 
be possible, and We c is apparently not the right 
form of Weber number to be used as a similarity 
parameter. 

Insufficient information was available to determine how 
the remaining parameters, <f), 0 and b, affect ice shape or 
how closely scale values of these parameters might 
need to be matched to their reference values. 
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a more important scaling parameter than one based on model size and air properties. Scale values of both the Reynolds and Weber 
numbers need to be in the range of 60 to 160 percent of the corresponding reference values. The effects of variations in other similarity 
parameters have yet to be established. 
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